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Abstract 

In this work, we present an implementation of quantum logic gates and algorithms in a three 
effective qubits system, represented by a (I = 7/2) NMR quadrupolar nuclei. To implement 
these protocols we have used the strong modulating pulses (SMP). The various stages of each 
implementation were verified by quantum state tomography (QST). It is presented here the results 
for the computational base states, Toffolli logic gates, and Deutsch-Jozsa and Grover algorithms. 
Also, we discuss the difficulties and advantages of implementing such protocols using the SMP 
technique in quadrupolar systems. 
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I. INTRODUCTION 



The characterization of the state of a quantum system is one of the most important steps 
in quantum physics and, in particular, in quantum computing (QC) [THS]. Nuclear magnetic 
resonance (NMR) has provided unique methods for demonstrating the realization of quantum 
logical operations and characterizing the quantum state of spin systems. The great sucess 
achieved so far by NMR QC is related to its ability for controlling the dynamics of the spins 
through radio frequency (RF) pulses [6l [7J. This great control also enables the creation 
of pulse sequences for the determination of elements of the nuclear spins density matrix, 
called the Quantum State Tomography process (QST). The first NMR QST method was 
developed by Chuang et al. [8] and optimized by Long et al. [9] for systems of heteronuclear 
coupled spins 1/2. It consists basically in applying a set of specially designed rotations on 
the different spins and reconstructing the density matrix from the resulting NMR spectra. 
This method was later adapted for homonuclear coupled spins 1/2 JTUJ [11] and also for 
quadrupolar spin 3/2 systems [121417] . where nonselective pulses were replaced by transition- 
selective RF pulses. The disadvantage of such method is the long time required to perform 
a selective pulse, producing considerable unwanted evolution of the system during the QST 
process. Besides, the mathematical formalism used to manipulate the transition operators 
[T8H20] are quite complicated and the resulting equations are not easily generalizable for 
nuclei with nuclear spin greater than 3/2. Thus, a new method for tomography [21] , based 
on nonselective pulses [22], was proposed and generalized to any spin system. 

Previous works, that used nonselective pulses to implement QST, explored a spin 3/2 
quadrupolar system (an effective two-qubits system) showing the great acurracy of the tech- 
nique [121 EI] also for relaxation studies [23TT26] and simulation of quantum systems [27] . 
However, when the spin quantum number is increased, some technical problems during the 
implementation may appear. For example, during each pulse of the QST process, several 
experimental errors, like oscillating field inhomogeneities, pulse imperfections, quadrupolar 
evolution, among others, may compromise the method performance. 

We report the QST adaptation for spin / = 7/2 in Section II. Testing procedures, results 
that confirm its efficiency and simulations that took experimental limitations into account 
are shown in Section III. In Section IV we show the experimental data obtained through 
the use of Strong Modulating Pulses (SMP) [2H] and QST for creating computational base 
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states, logic gates and quantum algorithms. 



II. DENSITY MATRIX TOMOGRAPHY METHOD 

In this section we will revise the QST method in details, with special emphasis on the 
general features that makes the method suitable for use in systems with arbitrary spins. 

A. Magnetization and the deviation density matrix 

In a NMR experiment, the observables are the components of the macroscopic nuclear 
magnetization, which are proportional to the average values of the nuclear spin operators 
(I a ), a = x,y,z: 

M a = (fi a ) octr(p/ a ) (1) 

where p is the density operator. In the high temperature approxiamtion, hu>L *C k B T, p can 
be written as: 

P « ^ + " A P ( 2 ) 
n being the number of qubits, e = hui/hBT ~ 10~ 5 is the ratio between the magnetic and 
thermal energies of the system, and Ap the traceless deviation matrix. In a quadrupolar 
spin system, the time evolution of the system magnetization, M (t), is obtained considering 
the following propagator 

U = e~— (3) 

in which 

H = H Z + H Q , (4) 

Hz and Hq are the Zeeman and the quadrupolar hamiltonian respectively. This hamiltonian 
is diagonal on the I 2 and I z eigenstates basis. Thus, because only the deviation matrix in 
Eq.pl) evolves under the influence of the propagator U, 

Ap(t) = UApUl (5) 

Where we considered Ap(0) = Ap (see Appendix |A| for details). 
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The expression for the evolution of the transverse magnetization can be obtained from 
the expected value of the I + operator. However, for a reason that will be clarified later, an 
a phase factor will be considered: 

M(t)=tr{UApU j I + }e ia (6) 

In Appendix [A] we show that Ap can be expanded in the polarization operator T^ m basis as: 

Ap = a ^e<^ m - m,) d l m > m (-0) T l>m , (7) 

l,m,m' 

Thus, inserting Eq.Q into Eq.([6]) and using the trace properties, we have: 

M{t)= fl^^^^'C,™ (-0) tr {T W C/t J + U) e ia . (8) 

l,m,m' 

This expression can be simplified using the properties of the polarization operators TJ jm 
[2H] (see Appendix [B]), 

Pl,m'] r)S = [Ti, m >] rjS 5 rtS+m >. (9) 

Using that 

[I + ] Sjr = [I + ] Sjr 8 s , r+h (10) 

we obtain, for the delta product 

fis,r+lSr,s+m' = ^m',-1? (H) 

that simplifies the magnetization to 

M (t) = J2 ai, m e^> +1) ^d l _ hm (-9) KK [I + ] s , [Ti-i] rtS . (12) 



l.m 



Where A s = e~ J T^*, and the X*X r = e~ l r *■ st provides the oscillation with Bohr frequencies. 

The term d l _ lm {—0) is the reduced Wigner function, which gives the dependence on the 
nutation angle 9. Hence, it is possible to maximize the magnetization by choosing 6 that 



make d l _ lm maximum. Using the properties (B5), (B2), and the hermiticiy of the density 
operator, we obtain 

a ltm = (-l)"X-m, (13) 
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it is possible to make the index change m — > —m (the summation of m runs over a symmet- 
rical range, |m| < I), so the final magnetization expression is reached: 



M (t) = J2 <me l{l - m) ^ +ia d^ m (-9) £ KK [I + ] s , r H 



Lm 



From the definitions: 



(14) 



fs,r = A* A 



(15) 



and 



we may write 



Ss,r = E ^e^-K^f )+-4 )m (-0) , 



l.m 



(16) 



(17) 



(18) 



B. Density matrix tomography 

With the magnetization written in terms of a sum involving the contributions of different 
orders of density matrix coherences, it is possible to design a coherence selection scheme, i.e. 
to write the NMR signal in terms of a single coherence order. Thus, the coherence selection 
was done with the temporal average of several signals with the form 



_ 1 p 

Sr,s ^ ^ S r ^ s ((fin; 



Or, 



n=0 



Inserting Eq.(17) in Eq.(19), we obtain 



(19) 



JSp-1 



r n=0 l,m 

With the angles <fi n and a n parametrized as 



(20) 
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n 7r 

<b n = 2n 1 — , 

Vn N p 2 

(m' - 1) 

a n = 2nn K J , (21) 

i\p 

the summation over n simplifies to 

£ e 2^(m>- m) = Np5mmi ^ Np > 1+m , (22) 



and Eq. (20) reduces to 



S r>r+1 (m') = £ < m ,< m , (-0) [A z ] r>r+1 . (23) 

It becomes evident now why the phase a was inserted in Eq. ^ to permit the emergence 
of the Kronecker delta. Notice that has to be set to the correct value during the signal. 
Since the operator TJ^o oc I is not accessible through NMR experiments, the latter equation 
will be rewritten with the auxiliar variable V: 
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r+1 



S r>r+1 (m') = al m ,d[ im , (-0) [Aj] rjH 

I' = {l,m'}. (24) 

This expression corresponds to the linear system 

AX = B (25) 

with 

[^Im-J'+i = j 

= <m'd[,rn> ("*) , (26) 

[S], = (m') , 

where 



6 



i = 1,2,.. .,25 



I = l',l' + 1,...,2S. 



(27) 



Therefore, for each average value, the respective coefficients a^ m i (with w! fixed) can be 
found. Because of the density operator hermiticity, only the coefficients for wl > need to 
be considered. To solve the system of equations, the least mean square method was used . 

III. EXPERIMENTAL TOMOGRAPHY TESTS 

In order to verify the efficiency of the tomography method for quadrupolar nuclei and 
the influence of experimental aspects in the method performance, some simple tests were 
performed. Those tests refer to the quadrupolar interaction influence of the sample, in- 
homogeneity of the radio-frequency field, and possible imprecision arising from the pulse 
calibration process. The sample used was cesium pentadecafluorooctanoate (CsPFO) and 
the nuclei observed was Cesium-133, wich has a spin of 7/2. The sample was dissolved in 
heavy water to the concentration of 37,5%. When the sample is placed in the magnet the 
liquid crystal miscels align with the field and, then, display an aerage internal electric field 
gradient. This will result in the following quadrupolar hamiltonian for the Cesium nuclei: 



The oo q is the quadrupolar frequency and its value is highly dependent on the sample 
temperature. Therefore, the sample temperature must be kept constant. The measurements 
were performed on a 400 MHz magnet using a 5 mm CP/MAS probe and a VARIAN Inova 
Unity spectrometer. 

A. Equilibrium, I y and even order coherence states creation 

In the base of Tj jTrv , the index m refers to contributions over the density operator for the 
m super-diagonal elements. For example, the m = 1 term has only the first superdiagonal 
terms non-null, independently of the rank /. The first test described is the acquisition of the 
equilibrium state I z , which has all off diagonal null elements (it is proportional to T 10 ). In 
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(3/, 2 - 1 2 ) ■ 



(28) 
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order to do so, we only need to keep the sample in the static magnetic field long enough for 
the spins to allign with it and then apply the tomography pulses. In the next test, the I y 
state is created by applying a (|)- x pulse on the equilibrium state. As I y have only the first 
(imaginary) sub and super-diagonal non-null elements, after the tomography sequence, only 
the spectrum corresponding to the first coherence (m = 1) will have non-zero intensities. 
The even order coherence test consists of creating states which have only pre-determined 
coherence orders by applying a pulse sequence composed of hard pulses and free evolution 
intervals. First, it is necessary to demonstrate that the free evolution changes the rank /, 
while the pulses change the coherence order m. To go from m to m', it is required that first 
the rank be changed from I to m', always obeying the constraints for the T^ m : 



m = —/,—/ + 1, ...,/ — 1, / 
I = 0,1, ...,25. 

1. Changing the rank 

For free evolutions, in the rotating frame on resonance, the evolution operator will be 

U e = e"* ,U Q = ^ [Ml - 1(1 + 1)1] . (29) 
Here, if the T^ m operator evolves under the effect of U e , the resulting T^ m will be given 

by 

f Um = U e T Wm Ul = Ki' T i',m. (30) 

V 

Thus, the information about the presence of the T;/ m component is in the b\y coefficient, 
wich may be evaluated by the projection of T^ m over T i m : 

Kl = tr {T^e-^T^e^} (31) 

For the rank transfer experiments, it is necessary to find the free evolution time t e for 
which the absolute value of b^L is maximum. For producing a state with only even order 
coherences this time is t P — — . 
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FIG. 1: Pulse sequence for the even order coherence creation 
2. Changing the order 

Using (— | )_ y pulses, based on the properties of rotations, and the fact that the density 



operator has contributions in the T/ m + Tj m format, we have (see Eq. (B9)): 



+ = [(-l) m d l , m (-e) + d l 0) _ m (-9)]T l , + 

+ E d U-m (-*) + (-1)"' d l m ,, m (-6)] (T hm , + Tl m ) , 

m'=l 

As showed by this expression, under unitary operations (i.e., rotations), the polarization 
operators T; m remains with the rank I unchanged. For the even order coherence state 
creation, it is necessary to find the theta angle value on the expression above for what the 
contributions by odd m are minimized (or null) and the even m maximized. By analysing 
the several contributions to the summation for m' = it was found that the even order 
maximizations occurs when 9 = | . In figure [T] the pulse sequence for the even order 
coherence creation is shown. The intermediate and final states are shown on figure |2j 
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FIG. 2: Real parts of intermediate £3 and final £4 states of the pulse sequence shown on figure [TJ 
All intensities normalized. 

B. Tomography Results 

As discussed before, the reduced Wigner functions d[ m (—9) are related to the nutation 
angle of each component and, in the tomography process, it is convenient that their values 
are as high as possible. Thus, the angle for which those functions are maximized should 



be the more convenient choice. However, even with the change to Eq. (14), that can lead 
to long pulses, making the pulse imperfections more severe and the tomography process 
less precise because selective excitation effects can start to show in the pulses mainly for 
large quadrupolar coupling. An alternative is to find angles that not necessarily maximize 
the d\ m (—0) function in particular, but those that for a given m make the values of these 
functions in relation to two or more values of I high enough for a considerable sensibility for 
this component. Thus, we used the tomography angles on the table [TJ which were selected 
based on this criteria. In the test analysis the following convention was adopted: FE = 
fidelity between the experimental state and its theoretical equivalent; FS = fidelity between 
the experimental state and the simulation. For the simulations, a e parameter related with 
the pulse precision was multiplied by the pulse intensity. A gaussian shape was adopted for 
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TABLE I: Nutation angle in radians for the tomography pulses 
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m=7 
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0.680 


1.020 
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1.094 
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0.268 




0.292 


0.426 


0.604 


1.404 
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0.730 


0.928 


1.426 



TABLE II: Experimental fidelities (in %) for each test. 





h 


ly 


Even order coherences 


FE 


99,66 


81,50 


75,25 


FS 


99,28 


97,73 


76,12 



the RF magnetic oscilating field, giving it a 5% maximum variation in intensity. 

The state with coherences of even orders was chosen as the guide for the simulations due 
to the fact that its expansions have components over a greater number of T^ m operators. 
With the FE determined, e and the quadrupolar frequency values were inserted on the 
simulation to reach a state with FS ~ FE. Therefore, the choice of these parameters was 
based on the experimental conditions, given a value of e — 0.95. When the parameters were 
defined, they were also used for the I z and I y simulations. 

The experimental results and the simulations are shown on the [3] and the fidelities are 
displayed on table [II} 

To acquire well-defined spectra, we chose for each coherence an average number obeying 
the progression 128, 144, 240. This gives the tomography procedure a duration on the 
order of 1 hour and 30 minutes. By a visual analysis of the experimental results and their 
respective simulations, it is possible to see that the matrix shape was well reproduced. 
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Experimental 



Simulation 




Real part Imaginary part Real part Imaginary part 



FIG. 3: At right, experimental data and, on left, simulation of the equilibrium state I z (i), Iy(H) 
and even coherence order state (iii). All intensities are normalized. 

IV. SMP RESULTS 



A. Computational base states using SMP 

Having confirmed the effectiveness of the tomography method, we now have a tool to 
determine completely the quantum state of the cesium nuclei in our sample. We wish 
then to demonstrate the possibility of performing a complete quantum algorithm in this 
system using the SMP technique [28]. However, before the implementation of the quantum 
algorithm we wish to demonstrate the construction of the pseudo-pure states used as the 
computational base. They will be the starting point of every algorithm. First we optimized 
the SMP pulse sequences for the preparation of the computational base states from 1 000) 
to | 111). The optmization was designed to employ four temporal averages, each one with 
ten segments. In Fig. [4], the experimental density matrices obtained through the quantum 
state tomography process are shown for each state, also with the corresponding simulated 
matrices that were expected. 



Table III gives the experimental fidelities corresponding to the density matrices presented 



in Fig. |4j The fidelities found during the optmization of the SMPs are very close to unity, 
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FIG. 4: Experimental density matrices and their respective expected simulated density matrices 

but as it can be seen, the experimental fidelities will always have a smaller value. That can 
be caused by experimental errors such as the ones described in previous section: a) pulse 
calibration; b) lack of spatial homogeneity of the B\ field of the RF pulses; c) limitations 
of the temporal, phase and amplitude resolution of the spectrometer; d) errors due to the 
response time of the RF coil during the rapid variations of power levels required by the 
SMP; e) relaxation processes during the pulse application. 

Among these different types of erros, there are some that are more influencing on the final 
result. From the simulations it was possible to determine which ones were more effective 
in reducing the fidelity of the measurements. The calibration errors are a strong factor in 
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TABLE III: Fidelity of the computational base states created 



State 


Experimental fidelity (%) 


1000) 


82,5 


1001) 


93,8 


|010) 


95,2 


|011) 


95,8 


1 100) 


92,8 


1 101) 


88,8 


|110) 


91,1 


1111) 


87,2 



the results. As demonstrated by the simulations, a 5% calibration error may result in a 
30% loss in fidelity. This shows the importance of a very precise pulse calibration. The 
inhomogeneity of the RF pulses is another important factor, but in this case not as much 
as the RF coil response time. The rapid variation of the SMPs amplitudes (sharp edges) is 
not sometimes well reproduced by the RF coil and this greatly reduces the fidelity. One way 
around this is forcing an SMP optimization that favors sequences with smooth amplitude 
transitions. 

An important thing to notice about the use of the SMP here is that the pulse durations 
for the creation of these states lie between 100 to 300 fis. If transition-selective pulses were 
used, the duration would be of the order of miliseconds. 

B. Logic Gates Implementation using SMP 

The next set of experiments we present consist of the logic gates implementation. For 
that purpose, we created a SMP that represents a logic gate and applied it to some of the 
computational base states. Then, we used the quantum state tomography method to check 
the output of the logic gate. The SMP pulse sequences for the gates were optimized with 25 
segments and the optmization process took about 10 minutes to reach fidelities over 99.99%. 
We chose to present here the results obtained for the Toffolli gates, since they are more 
graphically representative. 



14 



TABLE IV: Toffolli gate fidelities. 



Operation 


Fidelity (%) 


Toffollii |010) 


63,3 


Toffollii | 111) 


69,7 


Toffolli 2 |100) 


62,6 


Toffolli 2 1 101) 


66,2 


Toffolli 3 1 110) 


43,3 


Toffolli 3 1 101) 


60,3 



In Fig. [5] we show the tomographed density matrices before and after the application 
of the different Toffolli gates on the states |000), |001), |010), |100), |011) and |111). The 
experimental fidelities are shown in Table |IV| 




FIG. 5: Experimental density matrices of some computational base states before and after the 
application of the Toffolli gates with different target q-bits. 

As can be seen, despite the low fidelities, the final state can still be easily recognized as 
the one expected after the gate application. This demonstrates the correct performance of 
the SMP. To perform logic gates using transition-selective pulses on a / = 7/2 nuclear spin 
system, about 2 ms are required for the pulse duration. The SMP we used were at most 300 
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/xs long and therefore the 10 ms relaxation time for our sample has a negligible effect. 
C. Deutsch-Jozsa Algorithm 

With the use of SMPs it was possible to implement the Deutsch-Jozsa algorithm [30]. 
Here we present the result for two cases: one in which the tested function is balanced and 
one in which it is constant. In the case of a constant function the final result will be the 
state 1 000) or 1 001) and any other state if the function is balanced. In Fig. [6] the two results 
are shown. The constant function corresponded to one that always returned 1 for any input 
state. The balanced function returned for the states 1 00) and 1 10) and returned 1 for the 
states 1 1) e 1 1 1) . Thus, the final state of those constant and balanced functions is 1 001 ) 
and |011) respectively. We can visually recognize that the expected states are found in the 
experimental data, demonstrating the successful application of the SMP. 



Constant Function Balanced Function 




FIG. 6: Experimental density matrices obtained after the application of the Deutsch-Jozsa algo- 
rithm for a case of a constant and a case of balanced function. The final states should be 1 001 ) 
(constant) and 1 01 1) (balanced) 

Table [V] displays the experimental fidelity of the data obtained. The low fidelity values are 
attributed to the many experimental erros involved. Still, the final state is easily indentified. 
If conventional pulses were used, a sequence of at least 5 ms would be necessary to perform 
the algorithm, which would render it useless due to relaxation effects. Also a great analytical 
effort would be required to design a sequence that was the most efficient as possible. 
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TABLE V: Fidelity of the results obtained after the Deutsch-Jozsa SMP application 



State 


Experimental Fidelity (%) 


Deutsch-Jozsa f c 


64,3 


Deutsch-Jozsa fb 


65 



D. Grover Algorithm 

Our last result was the implementation of the Grover algorithm [3T]- This is an algorithm 
with practical uses, in which there is a yN order gain compared to the classical equivalent. 
Figure [7] shows the experimental data obtained for two cases of the algorithm. In the first 
case the search was for the |011) state and in the second case the search was for the 1 100) 
state. 



Search for 01 1 Search for 1 00 




FIG. 7: Experimental density matrices after the application of the Grover algorithm in search of 
the states |011) and 1 100) 

Theoretically, the final state of this algorithm for 3 qubits is not a computational base 
state. In theory, the final result has a 97% probability that it will return the state we were 
searching for. Amidst all our experimental errors, this is of little concern in our final results. 
On the experimental data, we can identify the states we were looking for even though the 



fidelities were not as high as desired. The fidelities are recorded in Table VI Due to all the 
experimental errors involved, these fidelities are as high as expected. 

We were able to check the performance of the SMP once again for another algorithm. 
Just as in the earlier cases, if we were to employ conventional pulses we would need a very 
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TABLE VI: Fidelities for the Grover algorithm 



State 


Experimental (%) 


Groveron 


52 


Groverioo 


52 



long pulse sequence (greater than the relaxation time) and also a great analytical effort to 
design it. 

V. CONCLUSIONS 

The QST method applied met the requirements of duration and precision of our system. 
That could be shown from the small differences between the simulated and experimental 
results. This could be evidenced by the tomography of the following states: 

• I z : The result is close to perfection, but this state has only coherences of order, and 
therefore it is not a good parameter for the tomography method analysis. 

• l y : Experimentally, it is possible to see the emerging of contributions on the real part 
with coherences of order 1 and, on the imaginary part, of coherences of order 2 (with 
much smaller intensities); yet, the most intense elements are exactly those that we 
expect theoretically. 

• Even order coherence: The elements that were theoretically expected can all be seen. 
The simulations with the introduction of the experimental imperfections, reproduced 
very well the real part observed experimentally. 

The construction of the computational base states was achieved with fidelity greater than 
80%, which is a good indication of the efficiency of the method. In the case of logic gates, 
the fidelities fall to the order of 60%, as expected due to the concatenation of SMP. In 
the same way, we find fidelities between 50% and 60% for the quantum algorithms. We 
have successfully shown the advantages and experimental applications of SMP for quantum 
information processing in 3 q-bit systems composed of 7/2 spin nuclei. The experimental 
data shows us that SMP allow us to perform logic operations in a time scale ten times smaller 
than with conventional pulses [To] . The implementation of complete quantum algorithms 
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was only made possible through the use of SMP. All the data was acquired through the 
quantum state tomography sequence which also had its efficiency tested. Both SMP and the 
quantum state tomography have many applications in NMR outside the scope of quantum 
information processing. SMP can make designing a pulse sequence a much simpler task 
without any analitical effort and sequences with much shorter duration and selectivity can 
be achieved through the use of SMP. 
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FIG. 8: Simplified description of the rotations 

Appendix A 



1. Rotations over the density operator 

To describe the rotations over the Ap operator, it will be expand in an ortonormal basis 
formed by the polarization operators {TJ )TO (S 1 )} [2U] : 

25 I 

Ap = ^2 ai , mTl > m - ( A1 ) 

1=1 m=—l 

Applying the rotation operator D (a, (3, 7) to the last expression, 

D(a,/3 n )ApD^(a,P,j) = ^ ai, m e~ imi ^ e~ im ' a d l m , m ((5) Ti jm /. (A2) 

l,m m' 

Any rotation of a coordenate system which takes (x, y, z) to (x', y', z') can be treated as 
a rotation of an angle Q around the n (0, $) axis, where G and $ are the same in both 
systems [29J , as illustred in the figure |8} 

The Euler angles (a, /3, 7) and the (Q, 0, $) angles are connected by: 



cos 



n 

2 

tan G 



cos 





2 

tan (f ) 



cos 



a + 7 



sm 







(A3) 
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Ti a — 7 

$ = — | -. 

2 2 

For a rotation axis constrained in the transversal plane, G = |. If the nutation axis 
makes a <fi angle with the x axis, $ = <ft. Then, for rotations of 9 about the axis Q = 8, we 
have: 



7 = 

a = - |, (A4) 
£ = -9. 

In NMR experiments, indicates the rotation associated to the pulse, the angle between 
the magnetic static field and the RF oscillating field (because they are perpendicular, 6 = |) 



and $ gives the phase of the RF field. Inserting Eq. (A4) in Eq. (B4) 



D l m ', m = e^-^-^d^i-e) (A5) 
Then, the expression of the operatorAp turned into A/5, is obtained: 

A ~P = E E V'^^^C.^ (-*) ZW (A6) 

I m,m' 

Appendix B 

1. Polarization operators T; m 

To tomograph the density matrix operator, it is necessary to expand the operator in 
a convenient basis. The polarization operators T^ m (S) were chosen for basis due to their 
properties under rotations. On the next calculus, all the algebric properties were extracted 
from [25] . The 7] >m operators are defined by the Clebsch-Gordan C l l '^ ni l2 m2 coefficients 



[Ti,m] cl>a - J 2S+1 C s',*,i,m> ( B1 ) 



where 
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\k-k\ < I <h + k, 

m = mi + m 2 = + 1, — 1, 1, 

I = 0,1,. ..,2S, 

a, a' = -S, -S + 1, S, 

a' = a + m, 



and 



(B2) 



The / index is the rank of the operator, where the m index is the coherence order. The 
rotations D (a, (3, 7) in the state space are descrited by the Wigner functions, D l m m ,. When 
the rotation operates in an eingenstate \l, m) of both J 2 and I z , it will result in a combination 
of states with a common value of I: 



D(a,p,-y)\l,m)= J2 D ™>,m l) % ™!) , 



m'=-l 



where a, (3 and 7 are the Euler angles and 



(B3) 



D l m ,, m (a,f3n) = e-« m ' a+m ^d l m , m (P) 
d l m , jm (p) = (km'le-^^m). 

The d l m , m ((5) are the reduced Wigner functions, which obeys: 



(B4) 



< m '(/3) = (-l) m ~ m '^ m ,- m '(/3)- (B5) 
The Clebsh-Gordan coefficients C l l '^ ni h m2 and the reduced Wigner functions d l m , m (/3) 
can be numerically calculated by the use of the following expressions: 



d l m >,m(P) = [(l + m)\(l-m)\(l + m')\(l-m')rx 

r / a\-\2l—2k+m—m' r . / a \ -i 2k— m+m' 

[cos (I) J L sm (f)J 



X 



k\ {I + m - k)\ (I - m' + k)\ (m' - m + k)\ ' 



(B6) 
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^aab 



X 



A (a, 6, C 

-i) 



E 



(a + a)! (a - a)! (c + 7)! (c - 7)! (2c + 1)' 
(& + /?)! (&-/?)! _ 

a - Q+2 (a + 6 - 7 - 2)! (6 + c - a - z)\ 



z\ (a — a — z)\ (c — 7 — z)! (a + b + c + 1 — z)\ ' 



x 



(B7) 



where 



A (a, 6, c) = 



(a + b- c)\(a-b + c)\ (-a + b + c)\ 



(a + b + c + l)\ 

Therefore, over a polarization operator, the rotation acts like: 



(B8) 



D (a, (3, 7) T l>m tf (a, 0, 7) = £ e^^d^m (/?) % 



(B9) 
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